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The  large  scale,  high  latitude  ionospheric  convection  pattern  is 
driven  by  the  interaction  of  the  solar  wind  and  the  interplane¬ 
tary  magnetic  field  with  the  earth’s  magnetosphere.  Ground-based 
digital  ionosondes,  with  spaced  receiving  antennas,  have  been  de¬ 
veloped  and  deployed  at  the  high  latitude  stations  of  Thule, 
Greenland  ana  Goose  Bay,  Labrador.  Results  from  these  stations 
demonstrate  that  the  ionospheric  convection  can  be  measured  for 
the  bottomside  ionosphere.  Data  from  Thule  consistently  show  the 
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predominant  antisunward  convection.  The  data  from  Goose  Bay  ind 
cate  the  sunward  return  flows  of  the  polar  plasma  convection  and 
the  switch  over  when  the  station  rotates  from  the  dusk  into  the 
dawn  cell.  These  data  also  illustrate  the  potential  for  system¬ 
atic  study  of  the  convection  patterns  that  is  possible  with  a 
network  of  ground-based  digital  ionosondes. 
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FOREWORD 


During  1986  two  international  conferences  were  held 
which  emphasized  solar  terrestrial  relationships: 

•  Solar  Terrestrial  Physics  Symposium 
XXVI  COSPAR,  Toulouse,  France 

30  June  -  12  July  1986 

•  International  Symposium  on  Space  Physics 
Beijing,  China 

10-14  November  1986 

A  number  of  sessions  in  each  conference  emphasized  develop¬ 
ments  in  understanding  of  the  polar  ionospheric  convection 
pattern  as  driven  by  the  interaction  of  the  solar  wind  and 
interplanetary  magnetic  field  with  the  earth's  magnetosphere. 

The  University  of  Lowell  Center  for  Atmospheric  Re¬ 
search,  with  cooperation  and  support  from  the  Air  Force  Geo¬ 
physics  Laboratory,  contributed  a  paper  at  each  conference. 
Both  papers  are  related  to  the  polar  ionospheric  convection 
and  the  systematic  study  that  is  now  feasible  with  a  network 
of  ground-based  digital  ionosondes. 

This  Scientific  Report  includes  the  papers  that  were 
presented  at  the  1986  international  conferences: 

Reinisch,  B.  W. ,  J.  Buchau  and  E.  J.  Weber,  "Digital 
Ionosonde  Observations  of  the  Polar  Cap  F  Region 
Convection,"  XXVI  COSPAR,  Toulouse,  France,  30  June 
-  12  July  1986  (to  be  published  in  Physica  Scripta, 
1987)  . 

Reinisch,  B.  W. ,  K.  Bibl  and  C.  G.  Dozois,  "Iono¬ 
spheric  Drift  Measurements  with  Digital  Ionosondes," 
International  Symposium  on  Space  Physics,  Beijing, 
China,  10-14  November  1986  (published  in  the  Sympos¬ 
ium  Proceedings,  1986)  . 
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DIGITAL  IONOSONDR  OBSERVATIONS  OF  THE  POLAR  CAP 
F  REGION  CONVECTION 

Bodo  W.  Reinisch 

University  of  Lowell  Canter  for  Ataospherlc  Raaaarch 

Jurgen  Buchau  and  Edward  J.  Weber 
Air  Force  Geophyslce  Laboratory 


ABSTRACT 


Ground-baaed  drift  obaervatlona  of  the  winter  polar  cap  F- 
reglon  ahow  that  the  aagnetospher Ically  Induced  lonoapherlc 
convection  can  be  aeasured  for  the  bottoaalde  lonoephere.  A 
digital  lonoaonde  with  four  epaced  receiving  antennae 
operated  at  Thule,  Greenland  (86°  CGL)  In  the  Doppler-drlf t 
aode .  A  nuaber  of  24-hour  aeasureaente  Indicate  that  the 
drift  direction  changee  linearly  aa  a  function  of  tlae  In 
accordance  with  the  predicted  antlaunward  convection  pattern. 
The  drift  velocltlee  vary  froa  300  to  900  a/a. 

Meaaureaenta  at  a  aubauroral  etatlon  (Gooae  Bay,  Labrador, 
63*  CGL)  with  the  eaae  epaced- antennae  -  Doppler  -  drift 
technique  ahow  a  steady  westward  drift  until  local  aagnetlc 
aldnlght  and  a  fast  switch-over  at  that  tlae  to  an  eastward 
drift.  We  conclude  that  the  observed  aubauroral  drifts  are 
the  eunward  return  flows  of  the  polar  plasaa  convection,  and 
the  switchover  occurs  when  the  station  rotates  froa  the  dusk 
cell  Into  the  dawn  cell. 


For  presentation  at  the  STP- International  Syaposlua  on 

Solar  Terrestrial  Physics 
XXVI  COSPAR,  Toulouse,  France 
i  30  June  -  12  July  1986 
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1.  INTRODUCTION 


Satellite  observations  have  established  the  exist¬ 
ence  of  a  two  cell  convection  pattern  in  the  polar  F  region 
drift  when  the  Interplanetary  magnetic  field  (IHF)  has  a 
southward  component  [1,  2].  This  convection  pattern  produces 
an  antlsunward  drift  at  the  highest  latitudes  and  sunward 
drift  at  lower  latitudes  in  the  auroral  zones.  The  plasma 
flow  pattern  is  less  predictable  and  more  irregular  when  the 
IHF  has  a  northward  component  [3,  4].  Uygant  et  al.  [S]  have 
shown  that  the  response  time  of  the  convection  to  changes  In 
the  IMF  is  two  hours  or  more. 

The  University  of  Lovell  and  the  Air  Force  Geo¬ 
physics  Laboratory  have  developed  ground-based  observation 
techniques  that  can  monitor  the  plasma  convection  as  a  func¬ 
tion  of  time.  tfe  have  begun  to  measure  the  F-region  drift  at 
three  high  latitude  stations:  (1)  Thule,  Greenland,  86°  CGL, 
(2)  Goose  Bay,  Labrador,  65°  CGL,  and  (3)  Argentia,  Newfound¬ 
land,  58°  CGL.  In  this  paper  we  -'ive  a  brief  description  of 
the  high-frequency  (HF)  radio  technique  used  to  measure  the 
ionospheric  drift  and  present  some  of  the  results  of  the 
Thule  and  Goose  Bay  observations . 


2.  SPACED -ANTENNA  DOPPLER  DRIFT  TECHNIQUE 

Spaced- antenna  HF  observations  of  ionospheric 
drifts  is  an  effective  method  to  study  the  dynamics  of  the 
ionized  atmosphere.  Vertically  transmitted  HF  waves  illu¬ 
minate  a  large  area  of  several  hundred  kilometers  diameter  in 
the  F  region;  an  array  of  antennas  receives  the  signals 
reflected  from  the  ionosphere.  For  the  meesurements 
described  in  this  paper  we  used  Digisondes  [6,  7],  i.e. 
edvanced  digital  ionosondes,  that  operated  alternately  in  the 
lonogram  and  the  drift  modes.  Actually,  the  ionograms  were 
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spaced  by  flva  ainuces  and  the  tine  In  between  was  filled 
with  a  number  of  18  sec  drift  measurements .  The  individual 
antennas  of  the  receiving  array  are  multiplexed  at  the  pulse 
repetition  rate  (200  Hz).  The  time  series  received  at  each 
of  the  antennas  is  Fourier  transformed  in  real  time  resulting 
in  four  complex  spectra  (in  the  case  of  four  receiving 
antennas)  at  the  end  of  each  measurement  period.  The 
spectral  resolution  was  0.125  Hz. 

Cross 'correlation  [8]  of  the  complex  spectra  from 
the  spaced  antennas  determines  the  angle  of  arrival  for  each 
spectral  component  containing  significant  signal  energy.  As 
a  result  of  this  analysis  one  can  construct  a  sky  map  showing 
the  location  of  each  reflection  point,  or  source  £  9 )  , 
specified  by  a  given  Doppler  frequency,  which  defines  the 
radial  velocity  component  of  the  moving  plasma  for  this 
source.  The  Doppler  frequency  d  is  given  by 

1 

d  -  -  j  v-k^  8-1,  2,  ...  (source  index) 

where  v  is  the  drift  velocity  and  k  the  wave  vector.  An 
example  of  such  a  sky  map  is  shown  in  Figure  1.  The  majority 
of  the  sources  are  located  in  the  south  east  with  a  near 
zenith  angle  of  about  30°.  The  Doppler  frequencies  vary  from 
♦  1.06  Hz  (labeled  1)  to  -0.81  Hz  (labeled  7).  By  assuming 
that  the  observed  Doppler  shifts  are  the  result  of  a  uniform 
bulk  motion  of  the  reflecting  plasma,  one  can  determine  the 
three*dlmenslonal  velocity  vector  v  which,  in  the  least 
squares  error  sense,  best  represents  the  Doppler  frequencies 
da  measured  at  the  source  locations  specified  by  kg . 

It  is  important  to  realize  that  the  14  sources 
shown  In  the  sky  map  of  Figure  1  existed  simultaneously 


within  a  n  18  second  tiae  window.  The  existence  of  the  aeny 
sources  was  first  established  by  the  spectral  analysis,  and 
then  their  location  was  found  by  cross -correlating  the 
antenna  signals  in  the  spectral  doaain.  The  diaensions  of 
the  receiving  antenna  array  are  shown  in  Figure  2  together 
with  the  array  pattern  for  a  10  MHz  signal.  It  is  evident 
that  the  angular  resolution  of  the  array  would  have  been  much 
too  low  to  resolve  the  different  sources  shown  in  Figure  1, 
if  it  had  not  been  for  the  preselection  by  the  Fourier 
analys is . 

The  cuaulative  sky  aap  in  Figure  3  shows  the 
strongest  sources  observed  during  11  drift  measurement 
periods,  identified  as  0,  1,  2,  ...  on  the  left  side  of  the 
figure.  While  the  reflection  points  (source  locations)  are 
changing  within  the  S  min  portrayed  here,  one  can  assume  that 
the  drift  velocity  vector  remains  constant  and  it  is 
therefore  possible  to  fit  a  velocity  vector  v  to  this 
accumulated  source  set.  It  is  easy  to  see  from  the  spatial 
distribution  of  the  Doppler  values  (shown  on  the  sky  map  on 
the  right  side  of  Figure  3)  that  the  drift  vector  at  this 
time  points  approximately  towards  the  west. 

3.  DRIFT  OBSERVATIONS  AT  GOOSE  BAY 

In  an  effort  to  study  the  diurnal  variation  of  the 
F  region  drift  we  calculated  the  average  drift  velocity  for 
every  IS  minute  interval.  The  direction  (azimuth)  and  the 
magnitude  of  the  horizontal  drift  component  are  plotted  as 
function  of  time  in  Figures  4,  5  and  6  for  three  winter  days 
at  Gooae  Bay.  In  all  three  cases,  the  drift  is  westward 
before  local  midnight  and  changes  within  about  one  hour  to  a 
predominantly  eastward  drift.  This  behavior  is  exactly  what 
is  expected  in  the  presence  of  the  two  cell  polar  convection 


pattern.  Gooaa  Bay  at  65°  CGL  la  undar  the  sunward  raturn 
flow  (aagnetlc  west)  of  tha  dusk  call  before  local  Midnight. 
At  aidnight,  Goose  Bay  rotates  through  the  Harang 
discontinuity  to  the  dawn  cell  where  the  sounder  now  senses 
the  sunward  flow  (aagnetic  eaat)  of  the  dawn  cell. 

The  typical  velocity  values  are  around  100  a/s,  but 
at  tiaas  the  velocity  goes  to  300  a/s  or  slightly  above  (not 
shown  in  the  figures).  The  vertical  velocity  coaponents  are 
indicated  in  the  speed  panels  by  +  (upward)  or  -  (downward) 
signs;  in  general  they  reaain  below  SO  a/s. 

Tha  exaaples  shown  here  were  selected  to 
deaonstrate  the  aidnight  velocity  reversal.  There  are  aany 
other  daya  where  the  reversal  occurs  but  is  not  as  clearly 
defined.  In  cases  where  the  velocities  are  very  saall  our 
aethod  produced  unreliable  results.  tfhat  needs  to  be  done  in 
the  future  is  to  increase  the  Doppler  resolution  by 
increasing  the  case  length  froa  18  sec  to  36  sec  or  aore. 

At  tiaes  the  drift  direction  changes  unexpectedly, 
as  at  07  AST  on  21  January  1982  (Figure  4);  for  about  three 
hours  the  drift  direction  is  westward  instead  of  eastward  as 
expected  in  the  aorning  hours.  Such  Irregular  changes  could 
be  the  result  of  an  expanding  and  shrinking  sice  of  the 
convection  pattern.  Vhen  the  pattern  expands  the  antlsunward 
flow  of  the  dawn  cell  aay  be  over  Goose  Bay  explaining  the 
observed  westward  drift. 


4.  DRIFT  OBSERVATIONS  AT  THULE 

In  winter  1983/1984  and  in  March  1985  AFGL's 
Airborne  Ionospheric  Observatory  was  deployed  at  Thule, 
Greenland  and  lonograas  and  drift  observations  were  conducted 
with  the  Dlglsonde  on-board  the  aircraft  and  four  crossed- 
loop  antennas  on  the  runway  nearby.  A  nuaber  of  24-hour  runs 
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are  shown  in  Figures  7,  8  and  9,  giving  examples  of  ths  thrss 
obsorvstionsl  periods.  Sines  Thuls  is  close  Co  che  center  of 
the  polar  cap  wa  expact  a  continuous  antisunward  flow  with 
its  direction  rotating  360°  in  24  hours.  In  magnetic 
coordinates  the  rotation  is  approximately  linear  in  time  as 
Thule  rotates  around  the  magnetic  pole.  In  all  the  examples 
shown,  the  general  trend  of  the  drift  direction  is 
approximately  antisunward.  On  9  December  1983  (Figure  7) 
there  are  some  large  fluctuations  in  the  direction  during  the 
evening  hours  and  they  coincide  with  the  occurrence  of  sun- 
aligned  arcs  observed  on  6300  A  with  an  all-sky  imaging 
photometer  [9].  Figure  8  shows  the  results  of  four  days  of 
observation  at  the  end  of  January  1984.  Again,  the  drift 
direction  closely  follows  the  antisunward  direction  which  is 
indicated  in  the  figure.  The  observed  velocities  vary  from 
about  300  to  800  m/s  which  is  three  times  higher  than  the 
typical  values  at  Goose  Bay.  Figure  9,  an  example  from  March 
1985,  again  showing  the  antisunward  plasma  flow. 


CONCLUSION 


The  exploratory  observations  made  at  a  polar  cap 
and  an  auroral  station  confirm  that  the  ground-based  spaced- 
antenna  Doppler  drift  technique  is  capable  of  monitoring  the 
high  latitude  convection  pattern.  The  average  plasma 
velocities  vary  from  100  to  300  m/s  at  the  auroral  station 
and  from  300  -  800  m/s  in  the  polar  cap.  Since  the  Danish 
Meteorological  Institute  is  now  operating  a  Diglsonde  at 
Qanaq  (87°  CGL)  and  the  Air  Force  Geophysics  Laboratory 
operates  Diglsonde  at  Goose  Bay  (65°  CGL)  and  Argentla  (58° 
CGL)  it  will  be  possible  to  more  closely  monitor  the  polar 
convection  and  to  correlate  the  observations  with  other 
geophysical  phenomena  (10,  11],  especially  the  geomagnetic 
activity,  and  with  the  variations  in  the  IMF. 
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IONOSPHERIC  DRIFT  MEASUREMENTS  WITH  DIGITAL  IONOSONDES 

Bodo  W.  Reinisch*,  Klaus  Bibl  and  Claude  G.  Dozois 
University  of  Lowell  Center  for  Atmospheric  Research 

ABSTRACT 

Motions  of  the  ionospheric  plasma  are  measured  by  Fourier  analyzing  the 
ionospheric  echoes  of  HF  transmissions.  The  transmit  antenna  of  a 
ground-based  ionosonde  uniformly  illuminates  the  ionosphere  in  a  cone  of 
about  30°  to  40°  around  zenith.  The  echoes  received  from  the  illumi¬ 
nated  ionosphere  have  a  Doppler  shift  u  -  2  k*v,  where  k  is  the  wave 
vector,  and  v  the  velocity  of  the  reflecting  plasma.  The  real  time  dis¬ 
crete  Fourier  transform  in  the  Digisonde  256  determines  w,  and  the  wave 
vector  k  (for  a  selected  sounding  frequency)  is  obtained  from  the  inci¬ 
dence  angle  of  the  echo  signal  with  Doppler  shift  w.  The  incidence 
angles  are  measured  with  a  receiving  array  of  seven  loop  antennas.  Dig¬ 
isonde  observations  at  high  latitude  show  very  consistent  drift  veloci¬ 
ties  reaching  magnitudes  of  900  m/sec  in  the  polar  cap  and  600  m/sec  in 
the  auroral  zone. 

INTRODUCTION 

In  the  past,  ground-based  ionospheric  drift  observations  employed  the 
similar  fading  (Briggs  et  al.,  1950;  Booker  et  al.,  1950)  or  correlation 
technique.  Here  we  discuss  the  more  direct  Doppler-drif t  technique 
(Pfister  and  Bibl,  1972)  developed  at  the  University  of  Lowell  for  use 
with  the  Digisonde,  a  modem  digital  ionosonde  (Reinisch,  1986;  Bibl  and 
Reinisch,  1978).  By  1988  more  than  thirty  Digisonde  256  stations  will 
be  operating  world-wide,  capable  of  providing  ionospheric  drift  informa¬ 
tion  for  the  E  and  F  regions.  It  seems  to  be  desirable,  therefore,  to 
describe  the  Digisonde  Doppler  Drift  (DDD)  technique  in  some  detail  and 
present  some  results  illustrating  the  potential  that  it  offers. 

DIGISONDE  DOPPLER  DRIFT  TECHNIQUE 

To  some  extent  the  DDD  technique  can  be  compared  with  the  well  known 
line-of-sight  Doppler  radar.  The  frequency  shift  of  the  reflected  radio 
signal  gives  the  radial  velocity  component  of  the  reflecting  target.  To 
measure  the  magnitude  and  direction  of  the  target  velocity  requires  Dop¬ 
pler  observations  by  three  spaced  radars.  In  case  of  the  ionospheric 
observations  we  are  not  dealing  with  a  point  target  but  with  a  large 
distributed  reflector.  The  incoherent  radar,  working  in  the  VHF  or  UHF 
band,  measures  three  (non-orthogonal)  plasma  drift  components  by  sequen¬ 
tially  pointing  the  radar  dish  into  three  different  directions  and  as¬ 
suming  spatial  and  temporal  uniformity  of  the  plasma  velocity. 

For  HF  observations  with  ionosondes  it  is  practically  impossible  to 
transmit  a  narrow  beam,  and  generally  the  transmit  antenna  illuminates 
the  ionosphere  in  a  cone  of  about  40°  around  zenith.  For  a  perfectly 
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smooth  and  spherically  symmetric  ionosphere  only  one  Fresnel  zone  in  the 
zenith  will  contribute  the  greatest  part  of  the  energy  returned  to  the 
sounder,  and  only  a  vertical  drift  component  can  be  determined.  In 
reality,  ionospheric  isodensity  surfaces  are  almost  never  smooth  and 
spherical,  and  echoes  are  returned  from  several  reflection  areas  at  the 
same  time.  Figure  1  shows  the  reflection  points  during  a  5 -minute 
observation  period  at  1315  local  time  at  Erie,  Colorado,  a  mid- latitude 


skymap  Colorado  ♦.«  mh< 

*  march  ms  mis  list 


Figure  1 

station.  There  is  the  small  spread  of  the  main  reflection  region  near 
the  zenith  point  which  is  usually  observed  during  quiet  daytime  condi¬ 
tions.  In  addition,  there  are  returns  with  zenith  angles  of  more  than 
10°.  In  this  skymap  each  reflection  point  is  specified  by  its  north- 
south  and  east-west  coordinates  and  by  its  measured  Doppler  frequency  (A 

-  +  0.05,  B  -  +  0.15  Hz .  1  -  -  0.05,  2  -  -  0.15  Hz,  ...).  The 

Doppler  shift  imposed  by  the  moving  ionosphere  is 

1 

d  --  k  •  y  s  -  1,  2,  ...  (source  index) 

S  7T  S 

where  y  is  the  drift  velocity  and  k  the  wave  vector. 

Actually,  such  skymaps  are  constructed  for  every  10  second  observation 
period,  each  showing  a  number  of  reflection  areas  or  sources  (Dozois, 
1983)  .  The  cumulative  skymap  in  Figure  1  is  the  superposition  of  about 
20  individual  skymaps. 

With  the  assumption  that  the  drift  velocity  field  is  uniform  over  the 
area  of  the  skymap  it  is  possible  to  calculate  a  three  dimensional  vec¬ 
tor  v  that  best  (in  a  least-squares  sense)  reproduces  the  measured  fre- 
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quency  shifts  at  all  the  observed  source  points.  The  vector  components 
vx,  Vy  and  vz  are  found  by  minimizing  the  error 

r  1  -|2 

«  -  X  d  -  -  (k  v  +  k  v  +  k  v  )  . 

“  s  it  xs  x  ys  y  zs  z'J 

The  components  of  the  wave  vector  k  and  kz  are  given  for  each 

source  by  its  location  relative  to  the  sounder. 

At  this  point  it  must  be  explained  how  the  Digisonde  measures  the  inci¬ 
dence  angle  (source  location)  of  one  signal  in  the  presence  of  many 
ocher  sources.  The  Digisonde  uses  an  array  of  seven  (sometimes  four) 
receiving  antennas  with  a  maximum  spacing  of  100  m.  The  3  dB  beamwidth 
for  a  3  MHz  signal  is  approximately  20°,  certainly  much  too  wide  to  sep¬ 
arate  the  sources  shown  in  Figure  1.  Indeed,  no  reasonable  receiving 
array  can  separate  in  terms  of  direction,  the  simultaneously  existing 
sources.  But  as  a  result  of  the  motion  of  the  ionosphere  these  sources 
will  have  different  Doppler  shifts  which  can  be  separated  by  Fourier 
analyzing  the  received  signals.  The  Digisonde  perfoms  in  real  time  a 
discrete  Fourier  transform  on  the  signal  for  each  of  the  seven  antennas. 
Cross -correlation  of  the  seven  complex  spectra  determines  the  angle  of 
arrival  for  each  spectral  component  containing  significant  signal 
energy,  resulting  in  the  sky  map  discussed  earlier.  The  sequence  of 
processes,  first  Fourier  analysis  then  incidence  angle  determination,  is 
important.  Inverting  the  sequence  would  lead  to  nonsensical  results 
since  the  phase  of  the  composite  signal  (originating  from  different 
sources)  has  little  physical  meaning.  For  example,  the  five  sources 
that  are  encircled  in  Figure  1  were  observed  within  one  10  second 
observation. 

Figures  2  and  3  show  the  results  of  drift  observations  in  Goose  Bay  (65° 
CGL)  and  Thule  (86°  CGL) ,  respectively.  For  the  auroral  station  at 
Goose  Bay  (Figure  2)  ,  the  direction  of  the  horizontal  drift  component 
changes  from  mainly  westward  before  local  midnight  to  mainly  eastward 
after  midnight.  The  drift  velocity  stays  between  50  and  300  m/s, 
occasionally  reaching  600  m/s.  The  drift  direction  in  the  polar  cap 
(Figure  3)  varies  linearly  with  time  and  the  velocities  stay  generally 
below  900  m/s. 
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